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Abstract 

We study the conditions under which the critical behavior of the three- 
dimensional mn-vector model does not belong to the spherically symmetrical 
universality class. In the calculations we rely on the field-theoretical renor- 
malization group approach in different regularization schemes adjusted by 
resummation and extended analysis of the series for renormalization-group 
functions which are known for the model in high orders of perturbation the- 
ory. The phase diagram of the three-dimensional mn-vector model is built 
marking out domains in the mn-plane where the model belongs to a given 
universality class. 

PACS numbers: 05.50.+q, 64.60.Ak 

According to the universality hypothesis pQ, asymptotic properties of the criti- 
cal behavior remain unchanged for different physical systems if these are described 
by the same global parameters. The field-theoretical renormalization group (RG) 
approach [2] naturally takes into account the global parameters and derives prop- 
erties of critical behavior from long distance properties of effective field theories. In 
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the present paper we study the long-distance properties of the d = 3-dimensional 
mn-vector model which is introduced by the following effective field-theoretical 
Hamiltonian [3]: 

"1 Tl 7~L / |j, \ 2 

Him] = I ^ E [l v0i 2 +^i0i 2 ]+5 E {\r\ 2 Y+ v i E \r\ 2 ) } • « 

J Z a=l a=\ \a=l / 

Here, <f) a = (<p a ' 1 , Q ' 2 , . . . , <f) a,rn ) is a tensor field of the dimension n and m along 
the first and the second indices; u and v are bare couplings; /i 2 , is a bare mass 
squared measuring the temperature distance to the critical point. 

Depending on the choice for the parameters m and n, the mn-vector model 
is known to describe phase transitions of various microscopic nature. The 
choice n = 1 comprises a bunch of systems that are characterized by an 0(m)- 
symmetric order parameter, while the limiting cases n — > and n — > oo correspond 
to these systems exposed to the quenched [I] and annealed [S] disorder respectively. 
The choice m = 1, arbitrary n corresponds to the cubic model A separate 
interest is provided by the case m = 2, n = 2 describing [Zj helical magnets 
and antiferromagnetic phase transitions in TbAu2, DyC2 as well as by the case 
m = 2, n = 3 describing antiferromagnetic phase transitions in TbD2, Nd. 

All of the mentioned cases of the mn-vector model were subjects of separate 
extensive studies (see e. g. Refs. jHJ 01101 an d references therein). They led to a 
consistent description of criticality in the 0(m) and cubic systems. In particular, 
the precise estimates of the critical exponents of the cubic and of the random 
Ising model were established both within high-order expansions of the massive 
and minimal subtraction field-theoretical RG schemes [HI EI] . On the contrary, the 
cases m = 2, n = 2, 3 remain to be controversial. Using general non-perturbative 
considerations it was shown ^l] that the theory {1} belongs to the 0(2) universality 
class. On the other hand the perturbative field-theoretical RG approach yielded 
mixed data, neither proving nor rejecting this result [HIE]- 

The studies infer that an intrinsic feature of the theory JIJ is an interplay 
between the 0(k) ("trivial") universality class (with k being dimensions m, or 
mn) and a new universality class. In this paper we address two problems that 
concern the crossovers in the mn-vector model and still attract attention. Firstly, 
we aim to obtain a map of universality classes of the theory (0) in the whole 
plane m > 0, n > 0. Such an analysis has been performed so far in the one- 
loop approximation [H]. We base the analysis on the high- loop expansions for 
the RG functions of the model (JTJ and its special cases; in order to refine the 
analysis we exploit Pade-Borel resummation [T3~] HI] of the (asymptotic) series 
under consideration. Secondly, we focus attention on the case m = 2, n = 2, 3 in 
order to explain why the highest orders of perturbation theory have not allowed so 
far to resolve what universality class is realized in the theory. We perform analysis 
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in different perturbative schemes and show that only certain of them give reliable 
answer. 

We analyze the theory applying the field-theoretical RG approach [2j within 
weak coupling expansion techniques. In the approach, a critical point corresponds 
to a reachable and stable fixed point (FP) of the RG transformation of a field 
theory. A FP {u*, v*} is determined as a simultaneous zero of the /3-functions 
describing the change of the renormalized couplings u and v under RG transfor- 
mations and being calculated as perturbative series in renormalized couplings. The 
equations for the FP read: 

(3 u {u*,v*) = u*<p(u*,v*) = 0, 

(2) 

Pv(u*,v*) = v*^{u*,v*) = 0, 

where we have explicitly shown that the structure of the (3-i unctions allows their 
factorization for the effective Hamiltonian (P). We make use of both the dimen- 
sional regularization with minimal subtraction [TH| and the fixed dimension renor- 
malization at zero external momenta and non-zero mass (massive) [TS| schemes. 
More precisely, we rely on the expansions for the /3-functions that are known at 
d = 3 with the accuracy of six loops in the massive scheme J7j and with five loop 
accuracy for the cases of 0(m)-vector and cubic models in the minimal subtraction 
scheme (Refs. [IHj and [THI correspondingly). 

Technically, the Eqs. (J2J) for the FP can be solved in two complementary ways. 
A perturbative solution is obtained by an expansion of the FP coordinates in a 
small parameter (e = 4 — d, with d being the space dimension of the model [201, m 
the minimal subtraction or massive schemes, or an auxiliary pseudo-e parameter 
[2T] in the massive scheme) around the Gaussian solution {u* = 0,v* = 0}. Such 
a way formally guarantees that the structure of solutions for the FPs remains the 
same after account of higher-order contributions once it has been established in 
the one-loop approximation. An alternative method (the 3d approach) consists in 
the solution of Eqs. (J2J) numerically ^01 122] at a given order of perturbation theory 
and provides less control on a loopwise upgrade. 

Within the perturbative approach, the conditions on m and n under which the 
critical behavior of the mn-vector model (0) belongs to a non-trivial universality 
class are known as Aharony conjecture and read j^j: 

n c < mn < m c n, n > 1. (3) 

Here, n c and m c stand for the marginal dimensions of the cubic model and of 
the random m-vector model. The conjecture is based on the one-loop stability 
analysis of four FP solutions compatible with the Eqs. (J2J). At d < 4, these 
are the Gaussian FP G {u* = 0, v* = 0}, the FPs P (mn) K = 0, v* ^ 0} 
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and Po(m) {u* 7^ 0,n* = 0} describing theories with one 4 coupling and thus 
corresponding to the 0{mn) and 0(m) universality classes, and, finally, the mixed 
FP M {u* ^ 0, v* ^ 0}. It is the stability of the FP M that is necessary for the 
appearance of a new non-trivial critical behavior. 

The 3d analysis of the theory (0) is obscured by our observation that at some 
choice of m and n more than four solutions for the FP are obtained. To convince 
ourselves that some of them are not a by-product of application of resummation 
procedures we propose to use the following argumentation. According to the basics 
of the RG theory, at the upper critical dimension d = 4 any 4 theory is governed 
by the Gaussian FP [2]. Therefore, any non-Gaussian solution at d — 4 is out of 
physical interest. If such a solution survives at any d < 4 and particularly at d = 3, 
we find natural to consider it physically meaningless by continuity. The situation 
becomes less clear if a FP cannot be continually traced back to a certain solution 
at d = 4 because it disappears at some 3 < d c < 4. In this case the stability of 
the estimate for d c as well as of the FP coordinates against application of different 
resummation procedures in different orders of perturbation theory might serve the 
purpose. It is to note here, that the special case of the theory (|T} with m = 2, n = 2 
is known to have exact mapping onto the model describing non-collinear magnetic 
ordering JB|. Within the massive RG scheme, standard six- loop 3d analysis of this 
model allowed to find a stable FP which does not have the counterpart within the 
perturbative e-expansion But one can not follow the evolution of the FP as 
d approaches 4 because in this case the resummation procedure is ill-defined j24l. 

We use both perturbative and 3d analysis as complementary ways to establish 
the map of universality classes of the theory ((TJ). We find that, in addition to 
the conditions (j3J), the high-order map is controlled by a degeneracy condition of 
one-loop equations for the FP |23] : 

n= 16(m-l) 
m(m + 8) 

Unlike order-dependent estimates for the marginal dimensions m c and n c , this 
equation is independent of the order of perturbation theory and is exact. We also 
observe that the results obtained with the account of high-order contributions differ 
qualitatively from those obtained in the one-loop approximation. We consider 
worth to mention three peculiarities. 

(i) We find a domain in the mn-plane where the high-loop resummed (3- 
functions produce no solution for the FP while such a solution exists in the one-loop 
approximation. In the mn-plane the domain spans from the vicinity of the point 
{m = m c ,n = n c /m c } upwards. There, we can solve the Eqs. (J2J) for the mixed 
FP reliably neither numerically at the fixed space dimension d = 3 nor by applica- 
tion of the pseudo-e expansion. In particular, though the pseudo-e expansion can 
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be formally constructed there, its analysis by means of Pade ^H] or Pade-Borel- 
Leroy ^3] technique produces highly chaotic values both for mixed FP coordinates 
and its stability exponents, (ii) We find a domain in the mn-plane where the 3d 
analysis reveals two solutions for the mixed FP M co-existing in opposite quad- 
rants of the -uf-plane. In the mn-plane the domain is located below the point 
{m = m c ,n = n c /m c }. Yet, we are always able to establish that one of the two 
solutions is unphysical in the sense explained above. The described phenomenon is 
quite stable with respect to the order of perturbation theory and to the type of the 
resummation procedure applied. In the perturbative approach only one solution 
for the mixed FP is present, (iii) We observe that a smooth change of parameters 
m, n in the mn-plane can show up as a complex abrupt trajectory of the FP M 
in the wv-plane. 

Realization of various universality classes of the theory (pQ) besides universal 
Eqs. ©-0 depends on non-universal initial conditions for couplings. Certain 
physical interpretations of the mn-vector model ((TJ impose restrictions for the 
signs of the couplings. Namely, a group including the cubic model (m = 1, Vn) 
and the cases m = 2, n = 2, 3 imply uq of any sign and vq > [HI Ej, whereas the 
microscopic base of the weakly diluted quenched m-vector model strictly defines 
u > 0, vq < Taking into account such a division along with the pseudo-e 
expansion based estimates for the marginal dimensions n c = 2.862 ± 0.005 9 
and m c = 1.912 ± 0.004 |2Zj, we arrive at the high-loop map of the universality 
classes of the theory Q as shown in the Fig. ^ There, the domains governed by 
different universality classes are bounded by lines for marginal dimensions and the 
degeneracy line. The FP M is stable for values m and n contained in dark regions. 
The stability regions of the FPs Po( m ) and Po(mn) ar e horisontally and vertically 
hatched. In the cross-hatched region in the Fig. both O(m) and O(mn) FPs 
are stable. Here, the choise of the universality class depends on the initial values 
of couplings u, v. They can be located in one of the two domains of uf-plane 
created by the separatrix, which is determined by the unstable mixed FP. The 
blank region in Fig. \I]p denotes the region of runaway solutions. Let us note that 
runaway solutions exist for the cubic-like models (Fig.^,), too; however, there still 
exist regions of initial couplings u, v starting from which the stable FP is attained. 

As we mentioned above the high-loop analysis of the theory encounters 
difficulties in some domains of the mn-plane. In particular, these are the domains 
where the mixed FP either disappears or can be given by two (physical and un- 
physical) solutions. Our direct calculations show that such domains (mainly) inset 
the regions where the FP M is expected to be unstable according to the Eqs. (JHJ)- 
(0J) and thus does not influence the analysis of the Fig. ^ However, even if the 
solution for the FP M is steadely recovered, its stability analysis is obscured for 
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Figure 1: The domains of FPs stability for the mn- vector model with different 
signs of couplings: Wu, v > (Fig. a) and u > 0, v < (Fig. b). The mixed FP 
is stable for values m and n from dark regions. The stability regions of 0(m) and 
0(mn) FPs are shaded by horizontal and vertical lines respectively. 



some values of m, n. In particular, the latter is observed for the physically inter- 
esting cases m = 2, n = 2, 3. At the rest of this paper, we aim to show that the 
reliability of the stability analysis depends on the choice of a series that is assumed 
as its basis. 

Indeed, the stability of a FP is governed by the condition 9?^) > with the 
stability exponents U{ being the eigenvalues of the matrix of derivatives = 
dfiui/ditj (ui = u,v) taken at the FP. For the case under consideration, m — 2, 
n = 2, 3, one of the eigenvalues (u^) is large and positive both at the FPs Po( m ) 
and M, so it is the sign of U\ that controls the stability of a FP. The exponent 
appears to be very small: an adjusted analysis of the 3d six-loop resummed RG 
expansions results in [S] Ui(m = 2, Vn) = 0.007(8) for the FP Po( m ) thus providing 
no definitive answer about its sign. The behavior of u)\ in different orders of 
perturbation theory can be explicitly demonstrated expanding the exponent at 
d = 3 in the pseudo-e expansion [21] parameter r up to the six-loop order: 

uj l (m = 2 1 \/n) = -l/5r+0.186074r 2 -0.000970r 3 +0.027858r 4 - 

0.014698r 5 +0.028096r 6 (5) 

and making an attempt to evaluate the exponent at r = 1 on the base of the Pade 
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table p]"]: 
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— u.zuuu 


— u.uioy 


n m aq 


u.uiou 


n nm 7 

— U.UUl ( 


u.uzu^ 


-0.1036 


-0.0149 


-0.0140 


0.0033 


0.0079 





-0.0717 


0.0251 


0.0053 


0.0209 








-0.0537 
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0.0113 











-0.0430 


0.0630 














-0.0351 


















In the table, numbers of the row and of the column correspond to the orders of 
denominator and numerator of appropriate Pade approximant for the exponent (J""}, 
the small numbers denote unreliable data, obtained on the base of pole-containing 
approximants, o means that the approximant can not be constructed. One can 
see, that the table shows no convergence even along the main diagonal and those 
parallel to it, where the Pade analysis is known to provide the best convergence of 
results [""""] . 

On the contrary, if one first defines a pseudo-e series for the value m = m c 
where the exponent o>i(m, Vn) changes its sign, one gets the series which has much 
better behavior [2*7] : 

m c = 4-8/3r + 0.766489r 2 - 0.293632r 3 + 0.193141r 4 - 0.192714r 5 . (6) 

Indeed, the corresponding Pade table for m c reads: 
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and leads to the conclusion m c < 2 already in the three-loop order (c.f. the 
convergence of the results along the diagonals of the table). A more efficient Pade- 
Borel-Leroy resummation procedure applied to the series (jBJ) results in an estimate 
[2*7] m c = 1.912 ±0.004. From here one concludes that iO\{rn = 2,Vn) > 0, the FP 
Po(m) at m = 2, n = 2, 3 is stable, and governs the critical behavior of the mn- 
vector model. In this way the perturbative RG scheme leads to the results which 
are in agreement with general considerations of Ref. "TT], where it was shown that 
the theory (jJQ) belongs to the 0(2) universality class for these field dimensions. 

Carrying out an analysis of conditions upon which the mn-vector model be- 
longs to the given universality class we met two problems which are worth to be 
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mentioned at the concluding part of this paper. The first is that an analysis of 
the resummed RG functions directly at fixed space dimensions may lead to an 
appearence of the unphysical FPs. One of the ways to check the reliabily of an 
analysis is to keep track of the evolution of the given FP with continuous change 
of d up to the upper critical dimension d — 4. The second observation concerns 
analysis of the FP stability: taking into considerations the contradictory results 
obtained by a direct analysis of the stability exponents we suggest that the most re- 
liable way to study the boundaries of universality classes in field-theoretical models 
with several couplings consists in an investigation of the expansions for marginal 
dimensions. We believe that these our observations might be useful at the analysis 
of critical propeties of other field-theoretical models of complicated symmetry. 

Work of Yu.H. and M.D. was supported in part by the Austrian Fonds zur 
Forderung der wissenschaftlichen Forschung, project No 16574-PHY and by the 
French-Ukrainian cooperation Dnipro project. 
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